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Introduction 

• Radiation transport codes, when combined with Risk 
Projection models, are main tool for shielding study 
and design. 

• Approaches to assess the accuracy of Transport Codes 

- Ground-based studies with defined beams and material 
layouts 

- Inter-comparison of transport code results for matched 
boundary conditions 

- Comparisons to flight measurements 

• NASA's HZETRN/QMSFRG code has a very high degree 
of congruence for each of these criteria. 





Approximate Composition 

^ 101 . 7 ^ 33 . 1^^36 

Density: 0.00194 g/cm 3 
Thickness: 1.2166 g/cm 2 
N: 2.09 10 22 atoms/g 
O: 6.81 10 21 atoms/g 
Al: 7.41 10 21 atoms/g 
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NSRL for Biophysics Applications 


Binary Filter 

0 ~ 24.225 g/cm 2 




Ion Chamber/SWIC (RW302) Booster Window 

0.09827 g/cm 2 0.10287 g/cm 2 



Air 

0.03482 g/cm 2 0.01471 g/cm 2 
0.08274 g/cm 2 


Air 

0.4257 g/cm 2 


Air 

0.12281 g/cm 2 


0.05379 ~ 0.08604 g/cm 2 


Heavy Ion Reactions 


Abrasion=projectile-target overlap 
(n, p, and cluster knock-out) 
Ablation = pre-fragment decay 
(n, p, d, t, h, alphas de-excitation) 
Coalescence = p and n knockout 
form bound states in couple phase 
space 
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Fragmentation Cross Sections: 
Comparison of QMSFRG to Si and Fe Beams 
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NSRL Bragg Curve Comparison to 
GCR Event-based Risk Model (GERM) 
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Poly Depth, g/cm 2 


Poly Depth, g/cm 2 


Cucinotta FA et al., Radiat Prot Dosimetry, 2011 



Thick Target Comparison with NASA's 
GERMCode* and GRNTRN Code* 


Iron (1 GeV/u) on Polyethylene 



10 V 

: Fe(5G,26) at 1053 MeV/amu on 10 g/cm^ Graphite -Epoxy 
6.0 mm Si detector 


Experiment (Normalized) 
NUCFRG2 
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*HZETRN uses identical Nuclear Cross Sections and Atomic Data 



Space Weather Prediction Center, NWS, NOAA 

NOAA/SWPC 12BIT 


Satellite Environment 


NOAA Scales Activity 
Range 1 (minor) to 5 (extreme) 

Past 24 hours 


NOAA Scale 


Current 


Geomagnetic Storms 
Solar Radiation Storms 
Radio Blackouts 


none 


none 


none 


none 


none 


none 
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Space Environmental Models 
Fit to Proton Measurements for Continuous Spectrum 

Functional Forms with Measurements 

• Exponential in Rigidity or Energy: ®(>R)=J 0 exp(-R/R 0 ) or ®(>E)=J 0 exp(-E/E 0 ) 

• Sum of Two Exponentials : ®(>E)=J 1 exp(-E/E 1 ) + J 2 exp(-E/E 2 ) 

• Weibull Function in Energy : ®(>E)=J 0 exp(- KE a ) 



Band Function with 4 Parameters (J^ y 17 y 2 , R 0 ): 
Double Power Law in Rigidity 


0(> R) = J 0 R~ n e~ RIR ° 




forR<(y 2 -y l )R 0 
h ~ n ^ ~ h ^ for R > (y 2 - y x )R 0 



E, MeV 


Flux, #/m**2 sr sec MeV/n 


Interplanetary Galactic Cosmic Ray Energy Spectra 

Advanced Composition Explorer/Cosmic Ray Isotope Spectrometer 
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Badhwar-O'Neill Model fit of ACE CRIS oxygen 
energy spectra measurements near solar 
minimum and near solar maximum 


Solar modulation parameter: 

ACE CRIS oxygen measurements (line); 
IMP-8 (Z>8) channel 7 measurements (o) 


O'Neil PM, 2010 



Geometry Models 

Shield Geometry Model and Shielding Analysis by CAD 


Structural Distribution Model for Layers of Spacecraft Using ProE™/Fishbowl 



Ray Tracing inside Spacecraft 


Color-coded Representation of Directional Shielding 






Human Geometry Models 
and Active Marrow Distributions 

Computerized Anatomical Male Male Adult voXel 
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Head and Neck 
12 . 2 % 

Chest 
26.1% 

Abdomen 
24.9% 

Pelvis 
33.4% 

Thighs/Upper Legs 
3.4% 


Lower Legs and Arms 
n/a 



MAlX 


Pelvic Region: 
20.9% 


H Flux (# particles/cm 2 MeV day) Neutron Flux (# particles/cm 2 MeV da 


Inter-Comparisons of Transport Codes 


n(C, AI20 g/cm 2 +H2O20 g/cm 2 ) 



10- 3 10- 2 10 -* 10 ° 10 1 10 2 10 3 10 4 10 s 10 6 
Energy (MeV) 




Energy (MeV) 



Energy (MeV) 


Energy (MeV) 

Heinbockel JH et al., NASA TP 2009-215560, 2009 




Comparisons with Flight Measurements 
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40 50 60 70 80 90 100 110 120 130 140 

Calculated dose rate. |i<3y day 



1.5-2.7X 
Albedo protons 
Albedo neutrons 
Secondary neutron 
Geomagnetic 
tansmission 
function 
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Calculated integral flux, number (cm sr day) 
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Badhwar GD, 1997 


Integral Flux, (cm 2 sr day) 



Evaluation of Detector Response 
- TEPC Response for Trapped Protons on STS-89 - 
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Phantom Torso Experiment (PTE) of ISS/STS 
TLD Dose Contours of Brain Slice 
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Badhwar GD et al., 2002 


Organ Dose Equivalent using CR-39/TLD, mSv 


Tissue 

Measured 

HZETRN/QMSFRG 

Difference (%) 

Skin 

4.5±0.05 

4.7 

4.4 

Thyroid 

4.0±0.21 

4.0 

0 

Bone surface 

5.210.22 

4.0 

-23.1 

Esophagus 

3.410.49 

3.7 

8.8 

Lung 

4.410.76 

3.8 

-13.6 

Stomach 

4.310.94 

3.6 

-16.3 

Liver 

4.010.51 

3.7 

-7.5 

Bone marrow 

3.410.40 

3.9 

14.7 

Colon 

3.610.42 

3.9 

8.3 

Bladder 

3.610.24 

3.5 

-2.8 

Gonad 

4.710.71 

3.9 

-17.0 

Chest 

4.510.11 

4.5 

0 

Remainder 

4.010.57 

4.0 

0 

Effective dose 

4.110.22 

3.9 

-4.9 


Yasuda et al., 2002 


Active Dosimetry Data, mGy/d 


Organ 

Trapped 

GCR 

Total 

Difference 


Expt 

Model 

Expt 

Model 

Expt 

Model 

(%) 

Brain 

0.051 

0.066 

0.076 

0.077 

0.127 

0.143 

13.3 

Thyroid 

0.062 

0.072 

0.074 

0.077 

0.136 

0.148 

9.4 

Heart 

0.054 

0.061 

0.075 

0.076 

0.129 

0.137 

6.7 

Stomach 

0.050 

0.057 

0.076 

0.077 

0.126 

0.133 

5.5 

Colon 

0.055 

0.056 

0.073 

0.076 

0.128 

0.131 

2.5 


Cucinotta FA et al., 2008 






E, mSv 



Predictions for Mars Mission 

1-y interplanetary space ■ 1-y Mars surface ■ 30-month Mars mission 
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Model-based Prediction of SPE Occurrence 
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SPE onset date 
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Model-based Prediction of SPE Fluence 


Propensity of SPEs: Hazard Function of Offset p Distribution Density Function 


m= 


/in 


+ 


K T(p + q) ( t 


4000 4000 T(p)T(q) 


V 


4000 


y-V t Y _1 

l 

J l 4000 J 


(0 <t< 4000 ) 


□ 95 percentile 
^ 25 percentile 


a 90 percentile x 75 percentile = Median 

x 10 percentile o 5 percentile The Carrington Event 
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Mission duration, d 


Effective dose on Mars Surface with MOLA Topography 


Altitude, 

km 

T,°C 

p, kPa 

Atmospheric shielding 
thickness, g/cm 2 
Low density High density 
model model 

8.0 

-41.16 

0.34 

0.14 

0.19 

4.0 

-34.99 

0.49 

6.73 

9.25 

2.0 

-33.00 

0.58 

10.97 

15.08 

0.0 

-31.00 

0.7 

16.00 

22.00 

-2.0 

-29.00 

0.84 

19.04 

26.17 

-4.0 

-27.01 

1.00 

22.64 

31.13 

-8.0 

-23.02 

1.44 

32.00 

44.00 


August 1972 SPE 







Conclusion 


Highly accurate descriptions of space environment models are available: 

> Inter-stellar GCR composition accuracy : ~5% for abundant elements (oxygen, carbon, 
and iron); less than 10% for all major GCR components; and solar modulation 
parameters with the 98.9% correlation in various spacecraft measurements. 

> Probabilistic SPE occurrence model as a tool for managing the risk. 

> Comprehensive catalogue of GLE fluences and spectra assembled for shielding 
design application using satellites and NM spectra; 

Radiation transport codes have been validated extensively: 

> QMSFRG model agrees for absorption a-section within +5% and elemental fragment 
a-section +25%. 

> Good agreement found from inter-comparisons of transport codes. 

> Comparison of model prediction to flight measurements: accuracy less than 15 % for 
GCR dose rates; ~25% for secondary particles ; and ±30% for quality factors by TEPC. 

> Minor scientific questions remained: low-energy light ion cross section, albedo 
protons, secondary pions, and kaons. 

Space Radiation Shield Design Tool for the reliable and realistic radiation simulation in the 
early design process of exploration missions: 

> Environmental models, shielding and body geometry models, atomic and nuclear 
interaction and fragmentation models are incorporated. 


